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We construct a density functional theory (DFT) for the sticky hard sphere (SHS) fluid which, 
like Rosenfeld's fundamental measure theory (FMT) for the hard sphere fluid [Phys. Rev. Lett. 
63, 980 (1989)], is based on a set of weighted densities and an exact result from scaled particle 
theory (SPT). It is demonstrated that the excess free energy density of the inhomogeneous SHS 
fluid $shs is uniquely defined when (a) it is solely a function of the weighted densities from Kierlik 
and Rosinberg's version of FMT [Phys. Rev. A 42, 3382 (1990)], (b) it satisfies the SPT differential 
equation, and (c) it yields any given direct correlation function (DCF) from the class of generalized 
Percus-Yevick closures introduced by Gazzillo and Giacometti [J. Chem. Phys. 120, 4742 (2004)]. 
The resulting DFT is shown to be in very good agreement with simulation data. In particular, this 
■ FMT yields the correct contact value of the density profiles with no adjustable parameters. Rather 

than requiring higher order DCFs, such as perturbative DFTs, our SHS FMT produces them. 
Interestingly, although equivalent to Kierlik and Rosinberg's FMT in the case of hard spheres, the 
set of weighted densities used for Rosenfeld's original FMT is insufficient for constructing a DFT 
which yields the SHS DCF. 
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I. INTRODUCTION 



, Density functional theory (DFT) was developed to characterize quantum systems in an external potential [l| and 
since the 1970's has been fruitfully brought to bear on classical fluids 0]. Indeed, here too it serves as an excellent 
O ! tool to study the microscopic structure of fluids under the influence of external potentials, such as confining walls, 
: and can serve as the sole prerequisite to map out complete phase diagrams that encompass various fluid and crys- 

| talline structures. Moreover, besides yielding a sound theoretical understanding of findings from colloidal physics and 
£SJ ■ computer simulations, DFT calculations can be extremely efficient computationally, especially if the systems under 
consideration allow for dimensional reduction, and hence in consequence largely outperform Monte Carlo (MC) sim- 
, ulations. In particular, DFT is especially convenient when phase diagrams or solvation free energies are calculated 
because it makes direct use of the grand potential f2. The basic formalism of DFT expresses f2 as a functional of the 
| | density profiles pi (r) of the different fluid components i = 1 , . . . , v as 

n[{ Pl (r)}} = J- id [{ Pl (r)}] + J- cx [{p,(r)}] + £ f dr Pi (r) [V^ (r) - W ] , (1) 

O »=i J 

' where fa is the chemical potential of species i, and I^ oxt (r) the external potential acting on particles of species i. The 
• , ideal gas part of the free energy is 

X: 

. S ; ^id[{Pi(r)}] = I dr Pl (r)(ln[A 3 p,(r)] - 1) , (2) 

i— 1 

where Aj is the thermal wavelength of species i and f3 = (fc^T) -1 . The excess part of the free energy •7 r cx[{pi(r)}] de- 
pends on the interactions between the particles and on temperature. For the given external potentials, the equilibrium 
density profiles are obtained by solving the system of coupled equations 8Cl/6pi = 0, i = 1, . . . , v. 

The above formalism is exact subject to the limitations of knowledge regarding J- ox for which we have only approx- 
imate descriptions for 3D fluids with non-trivial interactions. Among these approximations, Rosenfeld's fundamental 
measure theory (FMT) for the hard-sphere (HS) fluid Q provides an excellent framework that uses only exact re- 
sults for certain limiting cases which are generalized using dimensional analysis of a set of weighted densities. The 
theory applies equally to mixtures and, in part thanks to subsequent incorporations of highly accurate equations of 
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state [1, compares extremely well with computer simulations. The theory can be extended to treat crystallization 
of the HS fluid as well as the isotropic- nematic transition of non-spherical hard particles [9l4l^|. All of these 

developments make FMT a central and powerful tool for the study of entropy-governed hard-body fluids. 

Hard-body systems are an important test bed for theories which make use of a Weeks-Chandler-Andersen [l3| 
type decomposition of the interaction potential into a repulsive core and an attractive long-ranged tail. Moreover, 
they can often be considered as an accurate approximation to an entropy-dominated system. However, indeed many 
realistic fluids that occur in biological applications or colloidal science invariably exhibit attractive particle-particle 
interactions in addition to the hard-core repulsion. A simple model fluid with an attractive potential that has been 
extensively studied by theorists [3, [U| is the square- well (SW) fluid defined by an interaction potential 



a < r < a' (3) 




where a = 2R is the particle hard-core diameter, a' — a is the range of the attraction and esw is the well depth. 

The radial distribution function (RDF) g(r) and the direct correlation function (DCF) c(r) are related through the 
Ornstein-Zcrnikc (OZ) integral equation [16| 

g(r) - 1 = c(r) + p [ dr'c(r')(g(\r - r'|) - 1) . (4) 



However, the OZ equation can only be solved numerically within the Percus-Yevick (PY) closure for the potential 
0sw [13 ■ Alternatively, g(r) can be determined in computer simulations [l8j from which c(r) and thermodynamic 
properties of the SW fluid can be derived. However, these results cover only a limited selection of specific parameter 
choices for the model simulations and the results are clearly not in a convenient form for analysis. Hence, it is 
reasonable to seek a model with attractive interactions that would be solvable within the usual closure relations for 
the OZ equation. In the limiting case of surface adhesion where 

cr' - cr and (3e sw = In — — — oo , (5) 

12t((7 — a) 

Baxter's so-called sticky limit, with r _1 characterizing the strength of the attraction, allows the OZ equation to be 
solved within the PY closure fl9| . The results are convenient in a closed analytical form and compare well with 
computer simulations for various attraction strengths r _1 [2(J. Essential features of the sticky hard sphere (SHS) 
fluid, such as the existence of a critical point associated with a liquid-gas transition, are captured by the PY results 
[2l| . Subsequently, the model has been generalized for SHS mixtures where the PY closure leads to a set of coupled 
quadratic equations that can be solved analytically only in special cases [22|, [23[ . 

Far from being just a somewhat artificial theoretical construct, the SHS fluid has been demonstrated to provide an 
accurate description of colloidal systems with very short ranged attractions p4l - |26j and has been applied successfully 
to phenomena such as structural arrest in a micellar system (2~7j . competitive protein adsorption (28j . and even skim 
milk acidification [2!| . That said, the use of the SHS fluid to model experimental systems has been aimed at mapping 
out phase diagrams and calculating equations of state while the inhomogencous fluid has been studied less thoroughly. 
Firstly, this has to do with the difficulty of measuring structural properties such as density profiles and higher order 
correlation functions. Secondly, we lack a reliable tool to compute these quantities based on the analytical PY results. 
Existing DFTs for Baxter's SHS fluid are all based on an expansion of the one-particle DCF in powers of the deviation 
of the density profile from the bulk density. This expansion can only be performed systematically to linear order where 
it involves the DCF for which the PY results can be used [3(|. Moreover, the quality of the resulting DFT is rather 
poor, especially close to contact with a hard wall [3(j. Expanding the one-particle DCF to quadratic order requires a 
formula for the unknown three-particle DCF of the SHS fluid. The functional form of has been estimated by 
different authors based on assumptions regarding its range [3ll [32[ or its factorization in terms of c(r) [33]. In order 
to yield reasonable agreement with computer simulations all these extensions of the DFT by Choudhury and Ghosh 
[301 ] require an additional, density dependent, fitting parameter which is chosen such that the correct contact value of 
the density profile or the correct bulk pressure is guaranteed. 

Here we develop a DFT for the SHS fluid that is distinct from the above perturbative approaches. In the spirit of 
FMT we use a set of weighted densities that can be identified from a systematic inspection of the HS fluid in the low 
density limit. In order to construct a functional which yields the PY DCF for the SHS fluid, we supplement this with 
an exact relation from scaled particle theory (SPT) [34l . l35| . The main result is a unique construction for a broad 
class of closure schemes of which the PY closure is but one. No additional input and no fitting parameters arc used. 
Rather than requiring higher order DCFs our theory produces them. 
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In Scc.|TI]wc present the generalized PY closure as introduced by Gazzillo and Giacometti I36|. We then summarize 
in Sec. HIT] the FMT for the HS fluid in the version developed by Kierlik and Rosinberg (KR)j37|, HI] . In Sec.lPvHwe 
construct the new FMT for the single component SHS fluid, and then compare the results from the new functional 
with data from numerical simulations in Sec. |V] In the concluding section we place the new functional in a broader 
context, discussing both its merits and limitations with a special focus on future developments that might be pursued 
in the framework of our approach. 



II. GENERALIZED PERCUS-YEVICK APPROXIMATION 

Gazzillo and Giacometti introduced the generalized Percus-Yevick (GPY) approximation for the SHS fluid [3C 
The GPY is defined by a closure relation for the DCF c(r) of the form 



CGPY[r) 



-[l + 7(r)] 

^-shrink (t*) 3 



< r < a . 
a <r < a' , 



(6) 



0, 



r > & . 



Here 7(r) = pj dr'c(r')h(\r — r'|) is the convolution of the DCF with the total correlation function h(r) = g(r) — l, pis 
the particle number density of the bulk fluid. Depending on the choice of the function c s h r ink(0? when applied to the 
SHS fluid the GPY reduces to different familiar closure relations such as the mean spherical approximation (MSA), 
the modified mean spherical approximation (mMSA), and the PY closure [IB]. For example, choosing c s hrink(V) = 
/(r)[l +7(7")] in Eq. ©, where f(r) is the Mayer /-function, results in the PY closure relation cpyM = / (f)[l + 7( r )] 
for arbitrary r > 0. 

Gazzillo and Giacometti (3(| showed that in Baxter's sticky limit the GPY approximation allows for a general 
analytical solution of the OZ equation where the underlying closure relation enters solely through its approximation 
to the cavity function y a at contact. The resulting DCF is 



CGPY(r) 



r 



where 



1 



12r) 



(1 — i]) 2 I—?/ 
The PY closure corresponds to taking 

vFiv) = - 



-(a + b) 2 - aq a 



— - -a 2 \—J > Q(a - r) + q a crS(r - a) , 



b = - 



3rj 



2(l-77) 2 l-?/ 



and q n = - L ^— 
12 T 



(1-»7)t 



(7) 



(8) 



(9) 



(i-n)T 



3t* 



in which 



(i-v) 2 



(10) 



is the PY result for the HS fluid. The PY result for y a [Eq. ©] is exact up to linear order in 77, with a small deviation 
occurring at quadratic and higher orders. 

It turns out that the MSA leads to y^ SA = and hence does not capture the effect of surface adhesion in the 
SHS fluid and thus the results coincide with those for the HS fluid. The mMSA has y™ MSA = l which is the exact 
result in the low density limit; 77 — 0. Gazzillo and Giacometti (36j also introduce a CI approximation wherein 
7/^ :i = l + (| — - + y^pi) 7?, and hence it corresponds to the exact expression for y a truncated after the linear term 
in rj. 



III. FUNDAMENTAL MEASURE THEORY FOR THE HARD SPHERE FLUID 

Presently, FMT provides the most accurate free energy model for inhomogeneous HS fluids. Additionally, FMT 
^ives an excellent description of HS mixtures [39[ as well as fluids containing non-spherical hard particles f9j-fl2|] . The 



4 



theory originated with Roscnfcld's seminal paper Q and has been continuously modified by improving the underlying 
equation of state [1, [f| and the ability to capture the formation of HS crystals . There is an equivalent FMT 
formulation due to Kierlik and Rosinberg [37], [38j , which, in the case of the SHS fluid proves to be more useful than 
Rosenfeld's FMT (cf. Sec. IVI[) . In order to make this development reasonably self-contained, we summarize KR's 
version of FMT here. 

Fundamental measure theory provides an expression for the excess free energy .F cx which can be used in the 
framework of classical DFT to obtain the structural properties (density profiles and correlation functions) of fluids 
that are subject to external potentials, i.e., inhomogeneous fluids. As discussed above, we consider a fluid with v 
different components which have the corresponding density profiles Pi(r), i — 1, . . . , v. The FMT approach is based 
on weighted densities n a (r) written as 

n*(r)=ir / J dr' Pl (r'K (Q) (r-r') (11) 

i—l 

in terms of the weight functions u>\ a ^ of which there are four in KR's approach corresponding to four n a . These can 
be described in Fourier space: 

u4 3) (fc) = 47r[sin(fci? l ) - kR, cos^i?,)] A' 3 , (12) 

w f } (fc) = AirR, sm(kRi)/k , (13) 

a; 4 (1) (fc) = [sm(kRi) + kR, cos(kR l )]/2k , and (14) 

a4 0) (k) = cos(fci? l ) + (fcfli/2) sm(kR t ) , (15) 

or in real space as 

4 3 \r) = Q(Ri-r), (16) 

u J P(r)=6(R l -r), (17) 

u; 4 (1) (r) =d'(R l -r)/8n, and (18) 

4 0) (r) = -8"{Ri - r)/87r + 8' {R, - r)/2nR l - d{Ri - r)/2nR* . (19) 

Here, O is the Heaviside function and the primes denote differentiation. Using these weighted densities, the FMT 
excess free energy for the hard sphere mixture is readily obtained as 

= /dr$ HS ({na(r)}), (20) 

where the excess free energy density $hs is given by 

$ HS = -7i hi(l - na) + + ^ . (21) 

1 - n 3 247r(l - n 3 ) 2 

Fundamental measure theory, which may be derived independently of PY integral theory, can be shown to yield the 
same equation of state and the same DCF as obtained using the PY closure. 

IV. FUNDAMENTAL MEASURE THEORY FOR THE STICKY HARD SPHERE FLUID 

The FMT excess free energy density $hs obeys the SPT differential equation 

— =no -* + 5> (22) 

a 

which insures that the reversible work of adding a sphere to the fluid, in the limit of infinite sphere size, equals the 
volume of the sphere multiplied by the pressure of the bulk fluid [13, HH] . While this exact relation must hold for 
spheres with surface adhesion, there is no straightforward generalization of the principle for fluids with a finite range 
of interaction (in addition to the hard body repulsion) such as in the case of a SW fluid. For these fluids the question 
arises as to how the range of the interaction should scale as the size of the added sphere goes to infinity. However, 
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for the SHS fluid the surface attraction is subdominant with respect to the sphere volume in the large sphere limit. 
Therefore, Eq. (f2"2"j) must be obeyed. 

Requiring Eq. (|22|) to hold for the SHS fluid of radius R and Baxter parameter r implies that the corresponding 
free energy density $shs must have the form 

i , / Rn-2 \ R~ 2 n-) , ( Rn 2 \ , , 

*shs = *hs + -R _1 m 0i h — — - » 23 

V 1 - n 3 / 2?r \1- ra 3 / 

where 4>i and </>2 are arbitrary dimcnsionless functions of the dimcnsionless argument. We have made the additional 
assumption that $shs is linear in n\ which follows from inspection of the DCF cfmt obtained from an excess free 
energy density $ that is solely a function of the weighted densities n a . Then, in Fourier space we have 

CFMT(fc) = dfltr^^^^ ' (24) 

where, because we will treat only the single component fluid, we have dropped the second index of the weight functions. 
From the definition of the weight functions [Eqs. (fT2"|) - (fT5)) ] it can be seen that the inverse Fourier transform of 
a/ 1 )(/c)w( 1 ' l (/c) does not converge and hence, to insure that the DCF is defined in real space, we require $shs to be 
linear in ri\. A similar argument leads us to require $shs — ^hs to be independent of tiq. 

A crucial observation is that for any closure of the GPY class, there is a unique choice for the functions <fii and 
4>2 in Eq. |23|) such that cfmt [Eq. (|24|) using Eq. ([23]) ] equals cgpy from Eq. ([7]). This can be shown by comparing 
cfmt and cgpy in Fourier space. 

The functions 4>i and <f>2 are obtained through integration of the cavity function at contact y a associated with the 
given GPY closure. This is accomplished using the equations 

<j>' 1 {x) = -2y a {x)/x, (25) 

[^{x)]' = y^-xy a (x)+xy' a (x), (26) 

where x = and hence in the bulk fluid, where 71,3 = 7/ and n 2 = "irj/R, we have x = or ^ — The latter 
equation is used to rewrite y a = rjy a /T as a function of x rather than as a function of rj. The constants associated 
with the solution of Eqs. (|25[) and (|26[) are determined from the condition that 4>\ and 4*2 must vanish as x — > 0. 
The PY closure is achieved by using y^ Y from Eq. ^ in the formulas above and we find 

$ Y (x) =12t (VX- l) - Ax - 12rlnF + 2\/2(6t- l)lnZ, and (27) 



b\^{x) = (VX- 1) -6r 2 (5VX-3) +t(3 + 9x-Vx -xs/X 



x 

+ 36T 2 lnY+^(l-54T 2 )lnZ, (28) 
v2 

where 

Y = l(Vx + l)+^--^, and (30) 

(31) 

The MSA for the SHS fluid has y^ SA = 0, and hence <f>\ — 4>2 — 0, giving identical results to those for the HS fluid 
without surface adhesion. The mMSA and the CI approximation can be summarized by using the low density limit 
expression y^ DL = 1 + yii] with y\ = to recover mMSA and y\ = | — - + j^pj for the CI approximation. Solving 
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Eqs. ((25) and (J2BJ) with y^ DL yields the functions 

Vix 



(l + yi)ln (l + |) 



and (32) 



3 + a; 

(/>2 DL (x) t 6 + x + 6r ( 2 + x ) + 2 yi ( 9 + x + 3t ( 7 + 2a; )) + yf( 12 + x )] ln (* + 1 

- (2 + 4r + ri) - — j^i (18 + 42r + 5a; + 19™ + tx 2 ) 

Zt z Zt z {6 + X) 



III 



(216 + 126a; + I7x 2 ) . (33) 



12r 2 (3 + a;) 2 

The pressure jjshs related to $shs [Eq. ((23)) ] can be computed using either side of Eq. (f22j). viz., 

/?Pshs = «o - *shs + V = —z ■ (34) 

a 

With the advent of Eqs. (f2"5j) and ([2l>| it can be easily shown that pshs obeys the compressibility equation, which for 
the GPY closure is 



\d P J 



(35) 



T 



where a is given by the definition in Eq. © [for details see Ref. |36[. Therefore, the pressure underlying the FMT is 
always identical to the expression obtained from the GPY closure via the compressibility route. The quality of the 
mMSA, CI, and PY pressures is discussed in Ref. [H where the results are compared with computer simulations. 

V. COMPARISON WITH SIMULATION DATA 
A. Results for the radial distribution function 

As a first test for the consistency of the density functional f2, which uses the FMT excess free energy PJ- c -x[{Pi}] = 
J dr$sHs({"-a( r )}) with $shs from Eq. (|2"3"]l and the formulae 4>\ Y and 0f Y from Eqs. (|2T|) and (|28p. we minimize 
il in a spherical geometry. A center particle identical to the fluid particles is realized via an external potential and 
hence the resulting density profile corresponds to the RDF g(r) of the fluid up to a factor p. On the other hand, g{r) 
can be obtained directly from the DCF cpy [Eq. (0 with Eq. (J9j], using the OZ equation [Eq. ((4])]. Obviously, for 
an exact DCF and an exact functional the two results must be equal, but given that approximations are involved on 
both sides a deviation is to be expected. The degree of consistency and the agreement with numerical simulations 
help evaluate the quality of the FMT constructed here. 

For completeness, we include the formulae for the weighted densities in spherical geometry as they follow from 
Eqs. ([T2" )t to ([X5 ]) . The result is 



7T r +R 

n 3 {r) = - / drVp(r') [R 2 - (r - r') 2 } (36) 

r Jr-R 

n n rr+R 

n 2 (r) = — / drVp(r') (37) 

i rr+R r> 

"i ( r ) = ~r &r'r'p{r') + — [(r + R)p{r + R) + (r - R)p(r - R)} (38) 

4r J r _ R 4r 

Mr) = Y r t( r + R l 2 )P( r + R) + {r- R/ 2 )p(r - R)] 

-§[(r + R)p'(r +R)-(r- R)p'(r R)] , (39) 

where r, the distance from the center, is larger than R. Hence the problem of minimizing f2 in spherical geometry 
can be reduced to a computation in ID. 

In Fig.Q]we compare the FMT result for g(r) from the full minimization of Q to the results obtained (a) by using cpy 
in the OZ equation and (b) from MC simulations. Both the simulation data and the OZ results are taken from Ref. I20I. 
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FIG. 1: RDFs g(r) of SHS fluids with different Baxter parameters r and packing fractions r\. The full line represents the result 
from MC simulations [2(| and the dotted line was calculated in Ref. [2(J by using the PY result for the DCF, Eq. ((7) with 
Eq. (|9]), in the OZ integral equation. The dashed line is the result of the numerical minimization of the density functional fl 
which uses $shs from Eq. ([23) and the coefficients cf>\ Y and cj>l Y from Eqs. (|27j) and ([28j). a) r = 0.2 and rj = 0.32. b) r = 0.5 
and r] = 0.4. c) r = 1 and n = 0.5. 



For the three different choices of r and r\ the FMT and OZ results are consistent with each other and both agree 
well with the simulations. Deviations of the FMT results from the simulations occur mostly close to contact where 
the theory either underestimates (Fig. QJi) or overestimates (Figs. [T]d and c) the simulation data. Hence, there does 
not appear to be a systematic deviation and the good agreement of the two theoretical results with the simulations 
confirms the validity of the new functional. 



B. Results for the inhomogeneous fluid 



1. Density profiles 



A standard approach to test a DFT's ability to describe a fluid in an external potential which creates an inho- 
mogencity is to consider the fluid in contact with a planar hard wall. Due to the underlying symmetry the required 
numerical minimization of the functional f2 is effectively a ID computation and hence results are relatively easy to 
obtain. The formulas for the weighted densities in the planar case can be obtained by considering Eqs. (|36)) to ([39]) 
in the limit r — > 00. We compute density profiles for two different values of r and three different bulk densities pb 
using the FMTs based on the PY, mMSA, and CI closures for a fluid that is confined between two parallel walls at 
a distance of 1AR. In Fig. [2] the results are compared with data from canonical MC simulations taken from Ref. HO- 
The PY based FMT provides an excellent agreement with the simulations and the CI based FMT performs nearly as 
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FIG. 2: Density profiles of a SHS fluid confined between t wo p arallel hard walls at a distance 147?. The midplane is located at 
z — 0. Simulation data (open circles) are taken from Ref. WfJ. The curves are the results from the DFT presented in Sec. IIVI 
The different bulk densities pb of the fluid are given in the figure, a) Baxter's parameter is r = 0.2. b) Baxter's parameter is 
t = 0.5. 



well as the former. The mMSA based FMT shows deviations from the simulations which increase with decreasing r 
and increasing pb- Taking into account that (a) the mMSA FMT approximates the cavity function at contact y a (v) 
by J/(t(0) and (b) the CI FMT incorporates the density dependence of y a {ri) only to lowest order in r], we conclude 
that the quality of the underlying approximation for y&{r)) is not the crucial factor that determines the quality of 
the FMT. Indeed, the key ingredient appears to be the structure of the DCF cgpy(»") [Eq. ©] that follows from the 
GPY closure relation. The quality of the FMT results is clearly superior to the results obtained from the OZ equation 
based PY singlet approximation [4fl[4l| or the DFT by Choudhury and Ghosh I3C| . both of which deviate significantly 
from simulations especially close to contact with the wall. Higher order DFTs [3lTl33j require a fitting parameter in 
order to match the correct contact value or equation of state. In contrast, the fact that the FMT respects the contact 
theorem (42|, i.e. /3p — p c , without an adjustable parameter, insures that the density profile is nicely described by the 
present theory even close to the wall. 



2. Solvation force 

As a final test of how well the FMT compares with simulations we consider the solvation force / so i between two 
parallel hard plates at a distance d. The solvation force (per unit area) is closely related to the contact value p c of 
the density profile via f3f so \(d) = p c {d) — p c (oo). Here, we choose to plot f so i{d) rather than p c (d). This is because 
the accuracy of p c (oo) is assured by the underlying equation of state through the contact theorem (42[, while slight 
differences in the equation of state lead to different offsets for the p c (d) curves which obscures the comparison between 
the results. The solvation force for a fluid with Baxter parameter r = 0.2 and bulk density pb = 0.081 JR 3 is shown 
in Fig. [3^l. All the FMTs perform dramatically better than the OZ equation based singlet PY results [4fJ which are 
particularly poor for the contact density. As expected, the mMSA FMT is inferior to the PY and CI FMTs. In 
contrast to the density profiles in Fig. where the PY and CI FMT both capture the simulation data equally well, 
when describing / so i the PY FMT appears slightly superior to the CI FMT. This is consistent with the fact that a 
better approximation to y a underlies the PY FMT. Results for the more moderate case of a fluid which has r = 0.5 
and pb = 0.075/i? 3 are shown in Fig. The PY and CI FMTs show excellent agreement with the simulation data 
while the mMSA FMT are less accurate. The former versions of FMT clearly perform better than the DFT from 
Ref. [H (cf. figure 7 therein) . 



VI. CONCLUSION 

We have constructed a DFT for the SHS fluid which, in the spirit of Rosenfeld's FMT for the HS fluid Q , is based 
on a set of weighted densities which can be identified from a careful inspection of the dilute HS fluid. The main finding 
of this work is that using the FMT set of weighted densities (from the version of Kierlik and Rosinberg [13, [H|), in 
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d/R d/R 

FIG. 3: Solvation force / so i between two hard plates as obtained in simulations compared to the results from FMT. The 
solvation force is calculated from the contact values of the density profile according to /3f BO \(d) = p c (d) — p c (oo). a) The Baxter 
parameter is r = 0.2 and the bulk density is pb = 0.081/i? 3 . b) The Baxter parameter is r = 0.5 and the bulk density is 
p b = 0.075/i? 3 . 



combination with an exact relation from SPT [34|, l3fj , a density functional can be constructed which yields the PY 
DCF or the DCF resulting from any other closure of the broader GPY class [3(1 ■ The construction is uniquely 
determined by the set of weighted densities [Eq. (JTTJ with Eqs. - (fl9]l] the SPT differential equation [Eq. ([22]) ] 
and the given closure from the GPY class. 

We tested the performance of the new FMT functional for the SHS fluid by comparison with simulation data for 
the RDF and for density profiles of the confined fluid between two planar walls. The a gre ement is very good with the 
most significant improvement over previous theories such as the sing let PY results [H H or the DFT by Choudhury 
and Ghosh [3(| occurring in the density pro files close to contact with the wall. This reflects the property that the 
FMT complies with the contact theorem [42(, while other DFTs need an adjustable parameter to obtain the correct 
contact value [H, [Hj] . The results for the solvation force between two planar plates are also significantly better than 
the ones obtained with previous theories [HI, . 

Besides the obvious advantage of providing a very accurate description of the structure of the inhomogencous SHS 
fluid, the existence of an FMT formulation for the SHS fluid has several important implications for future developments. 
One aspect is the extension of the theory to SHS mixtures. Given that FMT provides a simple and accurate theory 
for HS mixtures, this new functional provides a starting point for the study of SHS mixtures. In particular, the 
weighted densities are the key to the simplicity of the HS mixture FMT. However, unlike this situation, the new 
functional depends on the particle radius R and the Baxter parameter r. Hence, a straightforward generalization of 
the functional with the help of weighted densities becomes impossible for mixtures in which R or r or both R and r 
are different. A future challenge resides in seeking an appropriate set of weighted densities for the SHS fluid which 
shifts the R and r dependence from the functional to the weighted densities, thereby opening the door for broad 
applicability to SHS mixtures. The complexion of such a set might be explored through incremental extension to 
simple binary mixtures some of which allow for an analytical solution within the PY closure [22l . |23| . Ultimately, one 
might be able to treat non-spherical sticky particles as was the case for hard particles Q . 

An additional future development includes extending the present DFT formulation to fluids with an attraction that 
has non-vanishing range such as the square-well fluid. Given that the PY closure cannot be solved analytically for 
these systems, it seems most likely that a successful road to an FMT would include the creation of new weighted 
densities which average the density profiles over a region that corresponds to the range of the attractive potential. In 
this context, studies of SW fluids with a very short but finite range of attraction [431 ] might provide a good starting 
point. 

Finally, there is the question whether the present FMT is suited for the study of SHS crystals. Problems such 
as spurious divergences for peaked density distributions for HS systems lead us to believe that the SHS FMT in its 
present form cannot properly describe a SHS crystal. However, these shortcomings have been successfully cured using 
different approaches within Rosenfeld's formulation of FMT Q. Nonetheless, these approaches have not yet 

been translated to the KR framework [37j which constitutes the basis of the SHS FMT. The alternative approach of 
formulating the present SHS FMT using Rosenfeld's set of weighted densities appears to be challenging. This follows 
from the observation that the 8 function in the DCF [Eq. ([7])] can only be generated using KR's weighted density 
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m. If the vectorial H2 in Roscnfcld's FMT, supplemented by 713 and 712, was used to construct $shs — ^hs the 
DCF in Fourier space would acquire a term proportional to w 2 (fc)uj2(fc) [see Ref. for the definition of w 2 ]- Such a 
formulation is unable to yield the term proportional to sinfci?cos kR/k originating from the 5-function in the DCF. 
This term can only be matched if the product LUi(k)ui2(k), with KR's ui(k) from Eq. (|T4"|). is used. Therefore, the 
extension of FMT to the SHS fluid provides an example where the equivalence of the KR and Rosenfeld formulations 
is no longer valid. This is an interesting and unexpected finding that illustrates the advantage of one formulation over 
the other depending on the task at hand. 
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